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Sr2ir04 has been suggested as a Mott insulator from a single J c g = 1/2 band, similar to the 
cuprates. However this picture is complicated by the measured large magnetic anisotropy and 
ferromagnetism. Based on a careful mapping to the J e g = l/2(pseudospin-l/2) space, we propose 
that the low energy electronic structure of Sr2lr04 can indeed be described by a SU(2) invariant 
pseudospin-1/2 Hubbard model very similar to that of the cuprates, but with a "twisted" coupling 
to external magnetic field (a g-tensor with a staggered antisymmetric component). This perspective 
naturally explains the magnetic properties of Sr2ir04. We also derive several simple facts based 
on this mapping and the known results about the Hubbard model and the cuprates, which may be 
tested in future experiments on S^IrCU. In particular we propose that (electron-) doping Sr2ir04 
can potentially realize high-temperature superconductivity. 

PACS numbers: 71.10.Fd,74.10.+v,75.30.Gw 



Various Ir oxides have recently become the platform to 
study the interplay between strong spin-orbit(SO) inter- 
action and strong correlation effects. There has been an 
experimental observation of a three-dimensional spin liq- 
uid in a hyper- kagome structure of Na4lr 3 0s [H . Theoret- 
ical proposals such as the realization of correlated topo- 
logical insulators [2] , the Kitaev model0|, and a Dirac 
scmimetal with surface "Fermi arcs" [J] in iridates have 
been made as well. Here we propose that doped Sr2lr04 
may realize high-temperature superconductivity similar 
to the cuprates. 

The crystal structure of Sr2lr04 consists of two- 
dimensional(2D) 1x02 layers, similar to the parent com- 
pound La2Cu04 of the cuprates. The main difference is 
that the oxygen octahedra surrounding Ir rotate along 
c-axis by about 11° in a staggered pattern, enlarging 
the unit cell by \/2 x \[2 x 2 [5|. The electronic struc- 
ture of Sr 2 Ir04 is quasi-2D, but is expected to have sev- 
eral differences from the cuprates. Ir 4+ has the elec- 
tronic structure 5d 5 , so the t2 g levels should to be ac- 
tive, while Cu 2+ with 3d 9 configuration has only the top 
e g level active. Ir as a bd transition metal is expected 
to have weaker correlation effects than 3d elements(e.g. 
Cu). At this point one may expect that Sr2lr04 is a 
(multi-band) weakly correlated metal. But strong spin- 
orbit coupling of Ir dramatically changes the story. The 
t2 g levels are split by SO interactions into a higher en- 
ergy Kramers doublet (the pseudospin-1/2 or J e g = 1/2 
states) and two pairs of lower energy ones[6|. These 
Jeff = 1/2 states are equal weight superpositions of all 
three t2 g orbitals, and this has been confirmed exper- 
imentally by resonant x-ray scattering^ and theoreti- 
cally by LDA+SO+U calculation Q. With d 5 configura- 
tion of Ir the J c ff = 1/2 states are half-filled. They have 
much smaller band width than expected for the tig lev- 
els without SO interaction and therefore have effectively 
enhanced correlation effect. In the end Sr2ix04 is a Mott 




FIG. 1: Schematic picture of one Ir02 layer. Large 
filled/open circles indicate the Ir atoms on two sublattices. 
Small open circles are oxygens. Small x, y are the global axis, 
while captial X, Y indicate local cubic axis(sublattices depen- 
dent). 



insulator and exhibits magnetic order below 24OK0- 11 1. 

It is then tempting to make the analogy between 
Sr2lr04 and the cuprates and speculate that doped 
Si - 2lr04 can also realize the interesting physics in doped 
cuprates, e.g. superconductivity, pseudogap, stripe for- 
mation, etc.. But strong SO interaction, different 
active orbitals and the rotation of oxygen octahedra 
seem to significantly complicate the problem. For ex- 
ample, Sr2lr04 has very anisotropic susceptibility and 
shows ferromagnetism(FM) with large ferromagnetic mo- 
ment ~ 0.14/ig per Ir[12j|. which was attributed to 
Dzyaloshinskii-Moriya(DM) interaction generated by the 
rotation of oxygen octahedra. However it has been 
pointed out by Jackeli and KhaliullinQ that the DM in- 
teraction can be removed by staggered rotation of pseu- 
dospin space on Ir sites. We will extend this considera- 
tion to the electronic model in a slightly different perspec- 
tive and discuss more details and consequences of this 
mapping. Finally we will argue that the analogy between 
Sr2ix04 and the cuprates can be established with care- 
ful interpretation, and that interesting physics of doped 
cuprates may also be realized in doped Sr 2 Ir04. 

The mapping to one band Hubbard model. To 
begin with we will treat Sr2lr04 as quasi-2D and consider 
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only one L-O2 layer, which is schematically illustrated in 
Fig. [TJ Label the rotation angle of oxygen octahedron 
around Ir site j by 9j = ej9, with ej = ±1 for the 
two sublattices and 9 s» 11° from experiments [5]. The 
crystal-field splitting of t 2g and e g levels and projection 
to Jeff = 1/2 states should be implemented in the rotated 
local cubic axis. Label the global axis by x, y, z and local 
cubic axis (on site j) by X, Y, Z (see Fig. H}. The unit 
vectors of these two coordinates systems are related by 

X — x cos6j+y sin#j, Y = —x sin9j+y cos9j, Z = z. 

(1) 

The J c ff = 1/2 states are (see e.g. Ref. [6j, the phase 
convention here is slightly different, i = 

|J c z ff = +1/2) = -L (+i\XYA) - \XZ,l)+i\YZ,l)) , 



I ^eff = -1/2) = ^= (-i\XY,i) + \XZ,t)+i\YZ,t)) ■ 



(2) 



XZ, YZ, XY are the t 2g orbitals defined in the local cu- 
bic axis, y, i indicate spin states (defined also in the 
local cubic axis). Note that although the elongation of 
oxygen octahedra along c-axis is expected to change the 
relative weights of the three orbitals [H this has not 
been observed in resonant x-ray scattering experiment 
or LDA+SO+U calculation @. 

As the first approximation, the effective electronic 
Hamiltonian should be the projection of full Hamilto- 
nian on the subspace of J c ff = 1/2 states. Considering 
first the Hamiltonian on the t 2g subspace, we expect the 
following, 1) the t 2g orbitals should be defined in the local 
cubic axis basis, because the crystal-field on Ir 5d orbitals 
from neighboring oxygens is diagonal only in the local 
cubic axis; 2) assuming that hoppings between Ir sites 
are mediated by the oxygen 2p orbitals, simple symme- 
try consideration shows that effective hoppings between 
nearest-neighbor Ir are orbital diagonal(one t2 g orbital 
does not hop to another orbital) only in the local cubic 
axis basis; 3) if the spin spaces are defined in the global 
axis basis, the effective hoppings of Ir t 2g orbitals will 
be real. Two tight-binding models on the t 2g subspace 
have been obtained by fitting LDA+SO+U dispersions 



in Ref. [8j and Ref. [13J, and both have this property 



of real orbital diagonal hoppings, but no clear interpre- 
tation was given. The discussion above shows that the 
orbitals in these models should be interpreted as the t 2g 
orbitals in the local cubic axis, while the spins in these 
models are defined in the global axis. 

The spin space on every site should be first rotated 
to local axis before the projection to the J c ff = 1/2 
states, because the spins used in ([2]) are defined in lo- 
cal axis. Namely we need to interpret the electron op- 
erators Cj a v used in these models, on site j for orbital 
a = XZ, YZ, XY with spin v, as creation operators for 



the states e ie " e j/ 2 |j, a, v) 
v =y, i respectively. 



where e v = ±1 for spin index 
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tions between t 2g orbitals will be projected into an onsite 
U term of the Hubbard model for the J e g = 1/2 states 
due to time-reversal symmetry and charge conservation. 

We take as a concrete example the tight-binding model 
of Ref. [13j. It involves nearest- neighbor(NN) XY hop- 
ping ti = 0.36eV, NN XZ(YZ) hopping along x{y) di- 
rection U = 0.37eV, NN XZ(YZ) hopping along y(x) 
direction t$ — 0.06eV, next-nearest-neighbor XY hop- 
ping t 2 = 0.18eV, and third-neighbor XY hopping £3 = 
0.09eV. The resulting one band Hubbard model after pro- 
jection is 

<jk>,a <<jk»,a 
«<jk»>,a j 

(3) 

with a =^, 1, and the effective hoppings are t = 
(1/3) (ii + U + t 5 )cos6» « 0.258eV, t = (1/3) (ii - 
U-t b )sm9 w -0.0045eV, t' = (l/3)t 2 ~ 0.06eV, 
t" = (l/3)t3 w 0.03eV. t is very small and will be ig- 
nored hereafter. In general t can be absorbed into t 
by a unitary transformation dj, a — ¥ e liea€j ^^ 2 dj ta with 
(j> — arctan(f/t), but we will not elaborate on this. The 
value of U has been estimated as ~ 2eV[l, |Tl|. This 
t — t' — t" — U model has been widely used as an ef- 
fective model for the cuprates, although the parameters 
here have different values. 

With large U and at half-filling the model © is an 
Mott insulator described by a pseudospin-1/2 model with 
SU(2) symmetry. If second- and third-neighbor t',t" 
are ignored the half-filling pseudospin model to the low- 
est order of t/U is just the Heisenberg AFM model of 
pseudospins J, H AFM = J2<jk> (^ t2 / u ) J j ' J fc- Eacn 
pseudospin has three components (a = 1,2,3) Jj t(l = 

(V 2 ) J2 a ,/3 d ],ai. (ja ) a i3 d j,i}^ wn ere a are Pauli matrices 
and a, (3 =t,| label the J* s = ±1/2 states. 

Coupling to external magnetic field. Although 
the effective model (0) looks exactly like the model of 
the cuprates, the coupling to external magnetic field in 
Sr2ir04 is quite different. 

Assume the coupling of magnetic field B on Ir 5d 
orbitals is described by the atomic form (more care- 
ful treatment can be found in, e.g., Ref. [6j), Hb = 
—fis^i ■ (L + 2S), where hb is the Bohr magneton. 
After projection to the J c ff = 1/2 states it becomes 
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H B = 2fi B B • J = 2hb(B x Ji + B Y J 2 + B Z J 3 ). Note 
that Bx,y,z are components of field on the local cubic 
axis, Bx = B • X etc.. Use the relation ([1}, the coupling 
on site j in terms of the field components on the global 
axis, B x ,B y ,B z , is 

+ Bj, y (Jj t 2 cos 6j + Jj.i sin 9j) + Bj >z Jj^.] 

(4) 

Therefore the observable magnetic moment Mj on site 
j has the following components on the global axis, 



M. 



j-y 



cosf -£, sin 








(5) 



If t in ([3]) is not ignored 9 sa 11° in ([5]) should be re- 
placed by 9 — arctan(f/i) ~ 12°. This nontrivial relation 
between moments M and pseudospins J, namely a ex- 
tensor with a staggered antisymmetric component, has 
several interesting consequences which we list below. 

• By quantum Monte Carlo studies [14| the square 
lattice Heisenberg model has a Neel ground state 
with staggered "magnetization" | ( e j J j ) | « 0.307. 
However because of the relation ([5]), the ordered 
moments do not form a simple collinear Neel pat- 
tern. If the ordered moments lie in the xy plane, 
they will be rotated together with the oxygen octa- 
hedra in a staggered pattern therefore create a net 
ferromagnetic moment per site, 2/i£-|(e :) Jj)|-sin# « 
Q.12^, B . This is very close to the experimentally 
observed value 0.14/zb [12 1. 



By the relation ([U we can relate pseudospin corre- 
lation functions of model ([3]) to moment correlation 
functions which is actually measured by suscepti- 
bility or magnetic neutron/x-ray scattering exper- 
iments. The Fourier components of moments with 
wavevector q and frequency u> is related to pseu- 
dospins by, 

Mq jUjX = -2^ B [cos6> Jq, w ,i - sin6» J q+ Q )Wl2 ], 

Mq,u,,y = -2(l B [cOs6 Jq, w ,2 + SU10 Jq+Q.w.l], 
Mq.u.z — — 2[l B Jq,uj,3- 

where Q = (tt,tt) is the wavevector of Neel or- 
der. In the paramagnetic phase the dynamical sus- 
ceptibility x ab (q,uj), which is proportional to the 
"moment structure factor" (M qiUia M^ q _ Ui (,), is re- 
lated to the dynamical pseudospin susceptibility 
X} b (q,o;) = <Ux/(q,w) oc (J q , w • J_ q _ w ) by 

= cos 2 9 x./(q, w) + sin 2 9 x,/(q + Q, w), 



X zz (q_,uj) = xj(q, w), and other components of x ab 
are zero. In particular the measured static uniform 
(oj = 0, q = 0) susceptibility in xy plane is actually 
a mixture of the uniform and staggered suscepti- 
bility of the SU(2) invariant Hubbard/Heisenberg 
model. This explains in a different perspective the 
measured large anisotropy of uniform susceptibil- 
ity and the ferromagnetic Curie- Weiss law 12|. In 
our picture the anisotropy is not mainly from easy 
axis interaction suggested by Ref. but from the 
mixing of large staggered susceptibility, and the 
FM Curie- Weiss law comes from the contribution 
of staggered susceptibility close to AFM Neel order 
of pseudospins. 

• In the high temperature paramagnetic phase above 
the Neel ordering temperature, the measured 
moment-moment correlation will be dominated by 
the staggered pseudospin correlation of a SU(2) 
invariant model, although the measured suscepti- 
bility shows significant anisotropy. The moment- 
moment correlation length will behave like the 2D 
Heisenberg model [15], which has recently been ob- 
served by magnetic x-ray scattering 16]. 

Possible high-temperature superconductivity. 

If the one band Hubbard model ([3]) is indeed a good 
approximation of the electronic struture of Sr2ir04, 
and if the high-temperature superconductivity in doped 
cuprates is indeed described by the one band Hubbard 
model, a natural consequence is that doped Sr2lr04 will 
realize high-temperature superconductivity. In the fol- 
lowing we list several direct consequences from this anal- 
ogy 

• It is believed that the sign and magnitude of t' is 
important for high-Tc in the cuprates and likely 
responsible for the particle-hole asymmetry of the 
phase diagram, (see e.g., Ref. 17]). The relative 
magnitude \t'/t\ w 0.23 for Sr2lr04 is similar to the 
cuprates. However the sign of t' for Sr2ir04 is op- 
posite to that of the cuprates. This can be remedied 
by a particle- hole transformation d- a — > e^dj a . 
Therefore we expect that the doping phase diagram 
of Sr2ir04 will be the particle-hole conjugate of 
the cuprates, in particular high-Tc will be easier to 
achieve on the electron-doped side of Sr2lr04, e.g. 
with La substitution of Sr. Interestingly electron- 
doped Sr2lr04_,5 has recently been synthesized and 
metallic behavior was reported for 8 = 0.04 [l8j . 

• The interlayer hopping of the cuprates is of the 
form = tj_QV 2 with v = (cosfej. — cosfc ?/ )/2, 
due to the d x 2_ y 2 orbital content[l9(. This to- 
gether with the d x 2_ y 2 nodal pairing symmetry sig- 
nificantly suppress transport along c-axis, making 
the superconducting properties of the cuprates very 
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anisotropy 
very small 
which 



anisotropic. However the resistivit" 
p c /p ab of Sr 2 Ir0 4 is only 10 2 - 10 3 [2( 
compared to 10 4 — 10 5 of the cuprates 21 
implies a larger t±Q for Sr 2 Ir04. The active or- 
bitals for Sr2lr04 is very different from the cuprates 
and the factor v 2 should be different and not van- 
ish on the nodal direction. Both facts suggest that 
Si - 2lr04 should have more isotropic superconduct- 
ing properties which is beneficial for practical ap- 
plications. 

• The pairing will be a pseudospin singlet d x 2_ y 2 
pairing and in many ways behave like the d-wave 
pairing of the cuprates. Phase sensitive and other 
indirect measurements used to determine the d- 
wave symmetry in the cuprates can be applied to 
doped Sr2ir04 as well. 

• The energy scale of the one band Hubbard model 
for Sr2lr04 is lower than that of the cuprates by 
about 50%. Therefore the T c of doped Sr 2 Ir0 4 will 
likely be lower than the cuprates. 

Discussion and Conclusion. The one band Hub- 
bard model ([3]) is of course the zeroth order approxima- 
tion of the low energy electronic structure of Si"2lr04. In 
real material the mixing between the J e s = 1/2 states 
and other states will generate anisotropy in pseudospin 
interactions, which will be important close to and below 
the Neel temperature. For magnetic properties above 
the Neel temperature and for electron-doped Sr2ir04 we 
believe this one band Hubbard model is still a good de- 
scription. 

The projection to one band Hubbard model was also 
implemented in Ref. The resulting hoppings re- 

ported in Equ. (7) (8) of Ref. [1] suggest that the au- 
thors of Ref. [8j interpreted the orbitals in their t2 g tight- 
binding model as the global axis basis xz,yz,xy. Here 
we have argued that the orbitals should be interpreted 
as the local axis basis which produces a projection result 
[i = -(2t /3) cos and h = 0] different from Ref. [|. 

In summary we have performed the projection of the 
electronic structure of Sr 2 Ir04 to the J c g = 1/2 states 
and carefully deduced the resulting one band Hubbard 
model and its interpretation. We provide another per- 
spective on the magnetic properties of Sr2lr04 by viewing 
it as a SU(2) invariant Hubbard/Heisenberg pseudospin- 
1/2 model, but with a twisted relation ([5|) between the 
observable moments and the pseudospin degrees of free- 
dom, namely a g-tensor with staggered antisymmetric 
component. One direct consequence is that the measured 
uniform susceptibility in ab plane is actually a mixture of 
uniform and staggered susceptibility of SU(2) invariant 
Hubbard/Heisenberg model. Despite the complication of 
strong SO interaction, different active orbitals and struc- 
ture distortion, the effective one band Hubbard model of 



Sr2lr04 remarkably resembles the cuprates. We thus pro- 
pose that doped Sr2ir04 can realize high-temperature su- 
perconductivity, and potentially other interesting physics 
of the cuprates. By comparing the model parameters 
we suggest that electron-doping of Sr 2 Ir04 will be the 
analogue of hole-doping of the cuprates. This can be 
achieved by La substitution of Sr, or O deficiency [lij . 
and maybe by field effect on thin films or interfacing 
with other oxides 0. We hope these simple theoretical 
observations will stimulate more experimental research 
on Sr2lr04. 
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